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We present an updated measurement of time-dependent CP asymmetries and the CP-odd fraction 
in the decay B° -> D* + D*~ using 232 x 10 6 BB pairs collected by the BABAR detector at the PEP-II 
B factory. We determine the CP-odd fraction to be 0.125 ± 0.044(stat) ± 0.007(syst). The time- 
dependent CP asymmetry parameters C+ and 5+ are determined to be 0.06±0.17(stat) ±0.03(syst) 
and — 0.75 ± 0.25(stat) ±0.03(syst), respectively. The Standard Model predicts these parameters to 
be and — sin2/3, respectively, in the absence of penguin amplitude contributions. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The time-dependent CP asymmetry measurement in 
B° — > D* + D*~ decay provides an important test 
of the Standard Model (SM). In the SM, CP vio- 
lation arises from a complex phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark- mixing matrix 0. 
Measurements of CP asymmetries by the BABAR and 
BELLE collaborations have firmly established this ef- 
fect in the B° — > J/ipK® decay [4| and related modes 
that are governed by the b — > ccs transition. The 
B° — > D* + D*~ decay is dominated by the b — * ccd 
transition. Within the framework of the SM, the CP 
asymmetry of B° — > D* + D*~ is related to sin2/3 when 
the correction due to penguin diagram contributions are 
neglected. The penguin-induced correction has been es- 
timated in models based on the factorization approxi- 
mation and heavy quark symmetry and was predicted 
to be about 2% |5j. A significant deviation of the mea- 
sured sin2/3 from the one observed in b — > ccs decays 
would be evidence for a new CP-violating interaction. 
The enhanced sensitivity of B° — > D* + D*~ to such a 
process arises from its much smaller SM amplitude com- 
pared with that of the b — > ccs transition. 

The B° -> D*+D*- decay proceeds through the CP- 
even S and D waves and through the CP-odd P wave. In 
this Letter, we present an improved measurement of the 
CP-odd fraction |fj, |j| R± based on a time-integrated 
one-dimensional angular analysis. We also present an 
improved measurement of the time-dependent CP asym- 
metry [fj 0] , obtained from a combined analysis of time- 
dependent flavor-tagged decays and the one-dimensional 
angular distribution of the decay products. 

The data used in this analysis comprise 232 million 
T(45) -> BB decays collected by the BABAR detec- 
tor at the PEP-II storage ring. The BABAR detector 
is described in detail elsewhere 0. We use a Monte 
Carlo (MC) simulation based on GEANT4 to vali- 
date the analysis procedure and to study the relevant 
backgrounds. 

We select B° — ► D* + D*~ decay by combining two 
charged D* candidates reconstructed in the modes 
D*+ -> D°ir+ and D*+ -> D+ir°. _ We in- 
clude the D* + D*~ combinations (D°ir + , D°n~) and 
(D°Tr+,D-ir°), but not (D+ir°, D--k°) because of the 
smaller branching fraction and larger backgrounds. To 



suppress the e + e~ — > qq (q = u,d,s, and c) continuum 
background, we require the ratio of the second and zeroth 
order Fox- Wolfram moments ^| to be less than 0.6. 

Candidates for D° and D + mesons are reconstructed 
in the modes D° — » K~it + , K~ir + ir°, K~tt + it + it~ , 
K° s ir+ir- and D+ -> K-tt+tt + , K° s ir+ , R-K+7T+. The 
reconstructed mass of the D° (D + ) candidate is required 
to be within 20 MeV/c 2 of its nominal mass except 
for the D° — » K ~ir + Tr Q candidate, where a looser require- 
ment of 40 MeV/c 2 is applied. 

The Kg candidates are reconstructed from two 
oppositely-charged tracks with an invariant mass within 
20 MeV/c 2 of the nominal K° mass. The \ 2 prob- 
ability of the 7r + 7r~ vertex fit must be greater than 
0.1%. Charged kaon candidates are required to be in- 
consistent with the pion hypothesis, as inferred from the 
Cherenkov angle measured by the Cherenkov detector 
and the ionization energyloss measured by the charged- 
particle tracking system [Sj ■ Neutral pion candidates are 
formed from two photons detected in the electromagnetic 
calorimeter 0, each with energy above 30 MeV. The 
mass of the pair must be within 30 MeV/c 2 of the nomi- 
nal 7T° mass, and their summed energy is required to be 
greater than 200 MeV. In addition, a mass-constrained 
fit is applied to the tt° candidates for further analysis. 

The D° and D + candidates are subject to a mass- 
constrained fit prior to the formation of the D* + can- 
didates. A slow 7r + from D* + decay is required to have 
a momentum in the T(4S I ) center-of-mass (CM) frame 
less than 450 MeV/c. A slow tt° from D* + must have a 
momentum between 70 and 450 McV/c in the CM frame. 
No requirement on the photon-energy sum is applied to 
the 7T° candidates from the D* + decays. 

For each B° — > D* + D*~ candidate, we construct a 
likelihood function |l2| £ maS s from the masses and mass 
uncertainties of the D and D* candidates. The likelihood 
£ maS s i s calculated as the product of the likelihoods for 
the D and D* candidates. The D mass resolution is 
modeled by a Gaussian whose variance is determined on 
a candidate-by-candidate basis. The D*-D mass differ- 
ence resolution is modeled by a double-Gaussian distri- 
bution whose parameters are determined from simulated 
events. The values of £ maS s and the difference of the 
B° candidate energy Eb from the beam energy -©Beam, 
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AE = E B - ^Bcam, in the Y(4S) CM frame are used 
to reduce the combinatoric background further. From 
the simulated events, the maximum allowed values of 

— ln£ mass and |A_E| are optimized for each individual 
final state to obtain the highest expected signal signif- 
icance using the previously measured B° — » D* + D*~ 
branching fraction |(|. 

The energy-substituted mass, toes = \/ ^Icam — P*b > 
where p* B is the B° candidate momentum in the T(45) 
CM frame, is used to extract the signal yield from the 
events satisfying the aforementioned selection. We se- 
lect the B° candidates that have toes > 5.23 GeV/c 2 . In 
cases where more than one B° candidate is reconstructed 
in an event, the candidate with the smallest value of 

— In £ maS s is chosen. A fit to the toes distribution with 
a probability density function (PDF) given by the sum 
of a Gaussian shape for the signal and an ARGUS 0] 
function for the background yields 391 ± 28(stat) signal 
events. In the region of toes > 5.27 GeV/c 2 , the signal 
purity is approximately 70%. 

In the transversity basis [l^. we define the following 
three angles: the angle 9\ between the momentum of the 
slow pion from the D*~ and the opposite direction of 
flight of the D* + in the D*~ rest frame; the polar angle 
0tr and azimuthal angle <fitr of the slow pion from the 
D* + defined in the D* + rest frame, where the opposite 
direction of flight of the D*~ is chosen as the x-axis, 
and the z-axis is defined as the normal to the D*~ decay 
plane. 

The time-dependent angular distribution of the decay 
products is given in Ref. [la] . Taking into account the de- 
tector angular acceptance efficiency and integrating over 
the decay time and the angles Q\ and tfitr, we obtain a 
one-dimensional differential decay rate: 

T a cos Otr 32n L 

X < — - — Io(cOS0 t r) H — /|| (cos 6* tr ) \ 

+ 2R ± cos 2 9 u . x /j.(oosfl t r)] , (1) 

where R±_ = \A^ 2 /(\A a \ 2 + \A\\ | 2 + |A ± | 2 ), a = (|A | 2 - 
l^||| 2 )/(l^o| 2 + |^||| 2 ), A Q is the amplitude for longitu- 
dinally polarized D* mesons, Am and Aj_ are the am- 
plitudes for parallel and perpendicular transversely po- 
larized D* mesons. The three efficiency moments, Ik 
(k = 0, ||, _L), arc defined as 

4(cOS0 tr )= fdcOsOi # t r Sfe(01,<£tr)£(01,0tr,<?f>tr), (2) 

where 50 = 4 cos 2 9\ cos 2 tr , 9\\ = 2 sin 2 9\ sin 2 (j> tv , g± = 
sin 2 61 , and e is the detector efficiency. The efficiency 
moments are parameterized as second-order even polyno- 
mials of cos#tr- Their parameter values are determined 
from the MC and are subsequently fixed in the likelihood 



fit to the differential decay distribution of cos#tr- In fact, 
the three Ik functions deviate only slightly from a con- 
stant, making the distribution, Eq.^ nearly independent 
of the amplitude ratio a. 

The CP-odd fraction R± is measured in a simultane- 
ous unbinned maximum likelihood fit to the cos 9t T and 
the Toes distribution. The background shape is mod- 
eled as an even second-order polynomial in cos^tr, while 
the signal PDF is given by Eq. ^ The finite detector 
resolution of the 9t r measurement is modeled as a dou- 
ble Gaussian plus a small tail component that accounts 
for misreconstructed events. The parameterization of 
the 9tr resolution function is fixed from the MC simu- 
lation and subsequently used to convolve the signal PDF 
in the maximum likelihood fit. Since the angle 9 tI is 
calculated with the slow pion from the D* + , we catego- 
rize events into three types: D* + D*~ — > (D°ir + , D°ir~), 
(D°Tr+,D-Tr°) 1 and {D+ir°, D°tt-), each with different 
signal- fraction parameters in the likelihood fit. Their 
angular efficiency moments and cos#tr resolutions are 
also separately determined from the MC simulation. The 
other parameters determined in the likelihood fit are the 
cos 9 tI background-shape parameter, three toes parame- 
ters (a and mean of the signal Gaussian, and the ARGUS 
shape parameter k), as well as R±. The fit to the data 
yields 

R x = 0.125 ± 0.044(stat) ± 0.007(syst). (3) 

The projections of the fitted result onto toes and cos #tr 
are shown in Fig. ^ 




5.23 5.24 5.25 5.26 5.27 5.28 5.29 5.3 -1 -0.8-0.6-0.4-0.2 -0 0.2 0.4 0.6 0.8 1 



m ES (GeV/c 2 ) cos9 |r 

FIG. 1: Measured distribution of m^s (left) and of cos#tr 
in the region m^s > 5.27 GeV/c 2 (right). The solid line is 
the projection of the fit result. The dotted line represents the 
background component. 



In the fit described above, the value of a is fixed to 
zero. We estimate the corresponding systematic uncer- 
tainty by varying its value from —1 to +1 and find neg- 
ligible change (less than 0.002) in the fitted value of R±. 
Other systematic uncertainties arise from the parame- 
terization of the angular resolution, the determination of 
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the efficiency moments, and the background parameteri- 
zation. The total systematic uncertainty on R± is 0.007, 
significantly smaller than the statistical error. 

We subsequently perform a combined analysis of the 
cos #tr distribution and the time dependence to extract 
the time-dependent CP asymmetry, using the event sam- 
ple described previously. We use information from the 
other B meson in the event to tag the initial flavor of the 
fully reconstructed B° — > D* + D*~ candidate. 

The decay rate /+(/_) for a neutral B meson accom- 
panied by a B°(B°) tag is given by 

/±(Ai,cos0 t r) oc e-l At l/ T * {G(l =F Aoj) T (1 - 2w) 

[Fsin(Am d At) + H cos (Am d At)]|, (4) 

where At = t ICC — t ta g is the difference between the proper 
decay time of the reconstructed B meson (B rcc ) and that 
of the tagging B meson (£>tag), tb° is the B° lifetime, and 
Arrid is the mass difference determined from the B a -B° 
oscillation frequency . The average mistag probability 
uj describes the effect of incorrect tags, and Aui is the 
difference between the mistag rate for B° and B°. The 
G, F and H coefficients arc defined as: 



G 
F 
H 



(1 - R±_) sin 2 tI + 2R± cos 2 9 tT , 

(1 - R±)S+ sin 2 tr - 2R ± S±_ cos 2 tr , (5) 

(1-R±)C+ sin 2 6» tr + 2R±C± cos 2 tr , 



where we allow the three transversity amplitudes to have 
different A fc = {q/p){A k /A k ) (k = 0, ||, _L) Q due to 
possibly different penguin-to-tree amplitude ratios, and 
define the CP asymmetry C k = 1 - | A fc | 2 / 1 + |A/c| 2 , S k = 
2S(Afc)/l + |Afe| 2 . Here we also have: 



C: = 



q||A||| 2 + C |A | 



\ A \\\ 2 



\Ao\ 



■j s+ — 



SAAu 



S \A Q \ 



\A 



(6) 



In the absence of penguin contributions, we expect that 
Co = C|| =C± = 0, and S = S\\=S± = - sin2/3. 

In Eq. the small angular acceptance effects are not 
incorporated, but absorbed into the "effective" value of 
R±, which is left free to vary in the final fit. No bias is 
seen in the resulting values of C + , Cl, S + , and S± in 
MC simulation. 

The technique used to measure the CP asymmetry is 
analogous to previous BABAR measurements as described 
in Ref . 0] . Only events with a At uncertainty less than 
2.5 ps and a measured |Af| less than 20 ps are accepted. 
We performed a simultaneous unbinned maximum like- 
lihood fit to the cos#t r , At, and mES distributions to 
extract the CP asymmetry. The signal PDF in 9 tr and 
At is given by Eq. The signal mistag probability is 
determined from a sample of neutral B decays to flavor 
eigenstates, -Bflav In the likelihood fit, the expression in 
Eq. |31 is convolved with an empirical At resolution func- 
tion determined from the -Bflav sample. The #tr resolution 
is accounted for in the same way as described previously. 



The background At distributions are parameterized 
with an empirical description that includes prompt and 
non-prompt components. We allow the non-prompt 
background to have two free parameters, C c ff and S e g, 
the effective CP asymmetries, in the likelihood fit. The 
background shape in tr is modeled as an even second- 
order polynomial in cos#tr, much as it is in the time- 
integrated angular analysis. 

The fit to the data yields 



C+ = 0.06±0.17(stat) ±0.03(syst), 

C± = -0.20±0.96(stat)±0.11(syst), 

S+ = -0.75±0.25(stat)±0.03(syst), 

S± = -1.75 ± 1.78(stat) ± 0.22(syst). 



(7) 



Fig. |21 shows the At distributions and asymmetries in 
yields between B° and B° tags, overlaid with the projec- 
tion of the likelihood fit result. Because the CP-odd frac- 
tion is small, we have rather large statistical uncertainties 
for the measured C± and Sj_ values. For comparison, we 
repeat the fit with the assumption that both CP-even 
and CP-odd states have the same CP asymmetry. We 
find that C+ = C ± = 0.03 ± 0.13(stat) ± 0.02(syst), and 
S+ = S± = -0.69 ± 0.23(stat) ± 0.03(syst). In both 
cases, the effective CP asymmetries in the background 
are found to be consistent with zero within the statisti- 
cal uncertainties. 




At(ps) 



FIG. 2: From top to bottom: the distribution of At in the 



region mES > 5.27 GeV/c 
the raw asymmetry (N B o - 



for B (B ) tag candidates, and 
Ng )/(N B o +Ng ), as functions 
of At. In the upper plot the solid (dashed) curves represent 
the fit projections in At for B° (B°) tags. 



The systematic uncertainties on C+, C±, S + and S± 
arise from the amount of possible backgrounds that tend 
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to peak under the signal and their CP asymmetry, the 
assumed parameterization of the At resolution function, 
the possible differences between the Bfl av and Bqp mistag 
fractions, knowledge of the event-by-event beam-spot po- 
sition, and the possible interference between the sup- 
pressed b — > tied amplitude and the favored b — > cud 
amplitude for some tag-side decays 0. It also includes 
the systematic uncertainties from the finite MC sample 
used to verify the fitting method. In general, all of the 
systematic uncertainties are found to be much smaller 
than the statistical uncertainties. 

In summary, we have reported measurements of the 
CP-odd fraction and time-dependent CP asymmetries for 
the decay B° — > D* + D*~ . The measurement supersedes 
the previous BABAR result , with more than 50 % re- 
duction in the statistical uncertainty, and indicates that 
B° — » D* + D*~ is mostly CP-even. The time-dependent 
asymmetries are found to be consistent with the SM pre- 
dictions within the statistical uncertainty. 
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